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Abstract 



We construct a four dimensional chiral M = 1 space-time supersymmetric 
Type I vacuum corresponding to a compactification on a toroidal Z2 ® Z2 
orbifold. Using recent results in four dimensional orientifolds, we argue that 
this model has a well defined world-sheet description. An interesting feature 
of this model is that the gauge group contains an SU(6) subgroup with three 
chiral generations. Moreover, this model contains D5-branes and therefore 
corresponds to a non-perturbative heterotic vacuum. This is the first ex- 
ample of a consistent chiral Af = 1 supersymmetric string vacuum which is 
non-perturbative from the heterotic viewpoint, has a perturbative description 
in a dual theory, and possesses some phenomenologically interesting charac- 
teristics. We also compute the tree-level superpotential in this theory. 
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One of the outstanding questions of string theory is whether it describes our universe. 
It is rather difficult to answer this question as the space of classical string vacua has a 
very large degeneracy, and perturbatively there lack objective criteria that would single 
out a particular string vacuum among the numerous possibilities. It seems reasonable to 
expect that non-perturbative string dynamics might lift, at least partially, this degeneracy 
of the moduli space. Understanding non-perturbative string dynamics is therefore of utmost 
importance. Moreover, in superstring theory supersymmetry is believed to be unbroken 
perturbatively. Since the world around us is not explicitly supersymmetric, understanding 
non-perturbative mechanisms of supersymmetry breaking in string theory is mandatory if 
string theory is to make contact with the real world. 

In recent years substantial progress has been made in understanding non-perturbative 
string dynamics. In particular, in ten dimensions there are five consistent superstring the- 
ories. The first four, Type IIA, Type IIB, E 8 ® E 8 heterotic and Spin(32)/Z 2 heterotic, 
are theories of oriented closed strings. Type I superstring is a theory of both unoriented 
closed plus open strings. Perturbatively, these five theories are distinct. Nonetheless, these 
theories exhibit a web of (conjectured) dualities which all seem to point to an underlying 
unified description. Most of these dualities are intrinsically non-perturbative and often allow 
to map non-perturbative phenomena in one theory to perturbative phenomena in another 
theory. 

Given the success of string dualities in shedding light on non-perturbative string dy- 
namics, it is natural to wonder whether one can find new consistent string vacua, which 
are inaccessible perturbatively, such that they possess at least some phenomenologically de- 
sirable qualitative features. Phenomenologically oriented model building in string theory 
for many years has been mostly confined to perturbative heterotic superstring framework 
where it was facilitated by the existence of relatively simple rules. In fact, many qualitative 
features such as realistic gauge groups, number of generations, etc., have been successfully 
obtained within this limited (as it is now believed to be) framework. It is natural to attempt 
construction of semi-realistic superstring vacua which would be non-perturbative from the 
heterotic viewpoint. Not only would this be encouraging as far as the hope for making 
contact with the real world, but it might also provide additional insight into the structure 
of string theory and, perhaps, give us clues how non-perturbative string theory can possibly 
solve some of the problems which might even seem hopeless from the perturbative heterotic 
viewpoint. 

In the light of the above discussion, one of the promising directions appears to be un- 
derstanding four dimensional Type I compactifications. One of the reasons to believe that 
this might indeed be a fruitful direction stems from the conjectured Type I-heterotic duality 
0. Thus, the tree-level relation between Type I and heterotic dilatons in D space-time 
dimensions || (which follows from Type I-heterotic duality in ten dimensions Jl|) reads: 

6-D , D — 2 . . . . 
<Ph = —^<Pi ~ -T7T log det 0, . (1) 
4 lo 

Here gj is the internal metric of the Type I compactification space, whereas (f>i and are the 
Type I and heterotic dilatons, respectively. From this one can see that (in four dimensions) 
there always exists a region in the moduli space where both Type I and heterotic string 
theories are weekly coupled, so one can hope to be able to understand the duality matching 
rather precisely. On the other hand, certain non-perturbative effects on the heterotic side 
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(such as dynamics of NS 5-branes which is difficult to study from the heterotic viewpoint) 
are mapped to perturbative Type I effects (e.g., D5-branes are the Type I duals of heterotic 
NS 5-branes) which are under much better control. 

The above expectations seem to be reasonable at the first sight, but there has been one 
circumstance preventing substantial progress in these directions. Thus, the four oriented 
closed string theories are relatively well understood as far as perturbative formulation is 
concerned. Conformal field theory and modular invariance serve as guiding principles for 
perturbative model building in closed string theories. Type I, however, lacks modular in- 
variance. Moreover, conformal field theories on world-sheets with boundaries (invariably 
present in open string theories) are still poorly understood. These have been some of the 
main reasons for lack of as deep understanding of perturbative Type I compactifications as 
in closed string theories. 

In the past years various unoriented closed plus open string vacua have been constructed 
using orientifold techniques. Type IIB orientifolds are generalized orbifolds that involve 
world-sheet parity reversal along with geometric symmetries of the theory Orientifold- 
ing procedure results in an unoriented closed string theory. Consistency then generically 
requires introducing open strings that can be viewed as starting and ending on D-branes 
H. In particular, Type I compactifications on toroidal orbifolds can be viewed as Type IIB 
orientifolds with a certain choice of the orientifold projection. Global Chan-Paton charges 
associated with D-branes manifest themselves as a gauge symmetry in space-time. The 
orientifold techniques have been successfully applied to the construction of six dimensional 
M = 1 space-time supersymmetric orientifolds of Type IIB compactified on orbifold limits 
of K3 (that is, toroidal orbifolds T 4 /Ztv, N = 2,3,4,6) 0||. These orientifold models 
generically contain more than one tensor multiplet and/or enhanced gauge symmetries from 
D5-branes in their massless spectra, and, therefore, describe six dimensional vacua which 
are non-perturbative from the heterotic viewpoint. 

The orientifold construction has subsequently been generalized to four dimensional M = 
1 space-time supersymmetric compactifications [flQ|,|2|j3~T^p~3]| . Several such orientifolds have 



been constructed. Some of these models, namely, the Z 3 [Efl, Z 7 [12| and Z 3 (8)Z 3 [ 1 3 1 orbifold 



models have perturbative heterotic duals 0-0]. Others, such as the Z 2 ®Z 2 [TO) and Z 6 [13 



orbifold models are non-perturbative from the heterotic viewpoint |ll|-Oj as they contain 
D5-branes. In particular, the Zq orbifold model of [OJ is the first known example of a 
consistent chiral M = 1 string vacuum in four dimensions that is non-perturbative from the 
heterotic viewpoint. 

Despite the above developments at some point it became clear that our understanding 
of orientifolds was incomplete. In particular, in some of the models discussed in [13, HI the 
tadpole cancellation conditions (derived using the perturbative orientifold approach, namely, 
via a straightforward generalization of the six dimensional tadpole cancellation conditions 
of Refs to four dimensions) allowed for no solutions. This, at least at the first sight, 
seems surprising as Type IIB compactifications on those orbifolds are well defined, and so 
should be the corresponding orientifolds. Clearly, a better understanding of the orientifold 
construction was called for. 



Recently, some progress has been made in this direction In particular, in [TJ 

conditions necessary for world-sheet consistency of six and four dimensional M = 1 su- 
persymmetric Type IIB orientifolds were studied. It was argued that in most cases orien- 
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tifolds contain sectors which are non-perturbative (i.e., these sectors have no world-sheet 
description). These sectors can be thought of as arising from D-branes wrapping various 
collapsed 2-cycles in the orbifold. In particular, such non-perturbative states are present 
in the "anomalous" models of Ref [14]] (as well as in other examples of this type recently 
discussed in Ref ||15|| ) . This resolves the corresponding "puzzles". Certain world-sheet con- 
sistency conditions in four dimensional cases (which are automatically satisfied in the six 
dimensions) were pointed out in [jTSJ] which indicate that the only four dimensional orien- 
tifolds that have perturbative description are those of Type IIB compactified on the Z 2 <8> Z 2 

Z 3 and Z 6 p"3| , and, finally, Z 2 <8> Z 2 <E> Z 3 (which we study in 



this paper) orbifolds. In particular, none of the other models considered in Refs [p| , |i5f have 
perturbative orientifold description, and even in the models with all tadpoles cancelled the 
massless spectra given in Refs [0,0 miss certain non-perturbative states. 

Although the number of perturbative Type IIB orientifolds in four dimensions appears to 
be rather constrained, it is important to understand all the cases at hand. The cases of most 
interest are those with D5-branes as they correspond to non-perturbative heterotic vacua. 
The only orientifold with perturbative description that has not been previously studied is 
that of Type IIB on the Z 2 <g> Z 2 ® Z3 orbifold. This is the case we will focus on in this paper. 
Happily, we will find a pleasant "surprise" , namely, the gauge group and matter content of 
this model (which is chiral) belong to the class of "semi-realistic" models which we were 
hoping to obtain to begin with. Moreover, this model is non-perturbative from the heterotic 
viewpoint (as it contains D5-branes), so it sheds some additional light on non-perturbative 
dynamics of heterotic NS 5-branes. 

The phenomenologically attractive features of this model are the gauge group and the 
number of chiral families. Thus, the gauge group of this model (for a particular configuration 
of D-branes) is [SU(6) <S> Sp(4)] 4 . The four copies of the SU(6) <g> Sp(A) gauge symmetry 
come from D9- and three sets of D5-branes. If we regard one of the SU(6) subgroups as the 
observable gauge group that contains the Standard Model, then the number of generations 
is three: there are 3 chiral 15's of SU(6) plus some number of 6's and 6's. (Note that 
the number of generations in SU(6) is determined by the number of chiral 15's.) Here we 
point out that in this model the charged matter is such that one cannot break the "grand 
unified" SU(6) gauge group down to the SU(3) C ® SU(2) W <g> U(l)y gauge group of the 
Standard Model by Higgsing. Nonetheless, it is very pleasing to observe that relatively 
simple Type I compactifications can yield certain qualitative features, such as the gauge 
group and the number of chiral generations, which are desirable phenomenologically. This 
is especially satisfying taking into account that this model is non-perturbative from the 
heterotic viewpoint. In fact, this is the first example of a chiral Af = 1 supersymmetric 
string vacuum which is non-perturbative from the heterotic viewpoint, has a perturbative 
description in a dual (that is, Type I) theory, and possesses some nice phenomenological 
characteristics. 

Clearly, better understanding of four dimensional compactifications of Type I on more 
generic Calabi-Yau spaces appears to be more than desirable. Construction of the model of 
this paper brings some hope that we might be able to learn more about phenomenological 
implications of non-perturbative heterotic vacua in the light of our discussion in the be- 
ginning of this section. Moreover, deepening or understanding of Type I compactifications 
would certainly be useful in understanding the underlying unifying theory and the nature 
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of string dualities. 



Type I on the Z 2 ® Z 2 <g> Z 3 Orbifold 

In the following we describe the construction of this Type I vacuum in detail. We 
compactify ten dimensional Type I superstring theory on a toroidal orbifold M. = T 6 /G 
where G is the orbifold group generated by three elements g,ri,r 2 . To describe the action 
of the orbifold group, let us for simplicity take T 6 to be a product of three two-tori: T 6 = 
T 2 ®T 2 ®T 2 . Let Zi (i — 1, 2, 3) be the complex coordinates parametrizing the T 2 's. Then 
the action of the orbifold group G is given by: 

gzi = ujZi , nzj = -(-l) SlJ Zj . (2) 

Here r 3 = rir 2 , and u = exp(27ri/3). Note that r\ generates a Z 2 group. Similarly, r 2 
generates a different Z 2 group. Finally, g generates a Z 3 group. Moreover, from @ it is clear 
that (7, n, r 2 commute, so that together they generate the orbifold group G = {g k , Tig k \i, k = 
1, 2, 3} ~ Z 3 ® Z 2 <g> Z 2 . (The element g 3 = 1 is the identity element of the group G.) 

It is convenient to view Type I compactification on M. as an orientifold of Type IIB 
compactified on M.. The orientifold projection is then simply given by Q, where Q is 
the world-sheet parity reversal that interchanges the left- and right-movers of Type IIB. 
Before we discuss this orientifold of Type IIB on A4, let us first understand Type IIB 
compactification on the M. orbifold. Note that M. is a Calabi-Yau three-fold with 577(3) 
holonomy. This implies that Type IIB on M. is a four dimensional theory with M = 2 
space-time supersymmetry. Let us compute the massless spectrum of this theory. For this 
it suffices to know the Hodge numbers of M. Calabi-Yau three-fold. Throughout this paper 
we will assume that the NS-NS antisymmetric tensor background is trivial (i.e., By = B-q = 
Bq = 0). Then it is not difficult to deduce the Hodge numbers in the orbifold sense. Here 
we give the contributions from each sector (untwisted and twisted) into the Hodge numbers 
(h 1 ' 1 , h 2 ' 1 ): untwisted: (3,0); g plus g~ l twisted (that is, Z 3 twisted): (9,0); twisted 
(that is, (Z 2 )j twisted): (6,0); gr, plus g~ x Ti twisted (that is, (Z 6 )j twisted): (2,0); total 
(h 1 ' 1 ,^' 1 ) = (36,0). 

Note that (h 1 ' 1 , h 2 ' 1 ) = (36,0) are the Hodge numbers of the Z-orbifold, that is, of the 
quotient M.' = T 6 /G', where G' = {l,g,g 2 } ~ Z 3 , and g acts as in (|2|) on the complex 
coordinates Z{ on T 6 . (Here and above we assume that the corresponding tori possess the 
symmetry by which we mod out or else the orbifold construction would be inconsistent.) 
Thus, M. and A4' topologically are the same. However, as we will see below, the perturbative 
orientifolds of Type IIB on these two quotients give rise to two different Type I models. (This 
is due to the different choices of perturbative (from the orientifold viewpoint) gauge bundles 
in the two cases.) 

Next, we turn to the Q orientifold of Type IIB on M.. Here some caution is required. 
The point is that the orientifold projection Q must be chosen to be the same as in the case of 
Type IIB on a smooth Calabi-Yau three-fold. Such an orientifold projection leads to Type 
I compactification on M.. The reasons why this choice of the orientifold projection is forced 
have been recently discussed at length in JL6|. In particular, we do not have an option of 



choosing the orientifold projection analogous to that in the six dimensional models of J7|||. 
On the other hand, the above Q orientifold projection is not a symmetry of Type IIB on 
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Ai at the orbifold conformal field theory point |16[ . The reason for this is that correctly 
reverses the world-sheet orientation of world-sheet bosonic and fermionic oscillators and left- 
and right-moving momenta, but fails to do the same with the twisted ground states. (Such 
a reversal would involve mapping the g and gri twisted ground states to the g~ x and g~ x Ti 



twisted ground states, respectively. In [ITH] such an orientation reversal was shown to be 



inconsistent.) This difficulty is circumvented by noting that the orientifold projection fl is 
consistent for smooth Calabi-Yau three-folds, and, in particular, for a blown-up version of the 
Ai orbifold. Thus, once the appropriate blow-ups are performed, the orientifold procedure 
is well defined. (In particular, the Z 3 and (Z 6 )j twisted sector singularities must be blown 



up. The (Z 2 )i singularities are harmless, however ||16|| .) 

Another important issue is related to the question of whether all the sectors of the 
orientifold have a well defined world-sheet interpretation. Note that the orientifold group is 
given by O = {g k , r i g k 1 flg k , Qrig k \i, k = 1, 2, 3}. The sectors labeled by g k , correspond 
to the (unoriented) closed string sectors. The presence of the element Q G O implies that 
we must introduce 32 D9-branes (just at in ten dimensional Type I). Similarly, the presence 
of the elements flri G O requires introduction of three sets of D5-branes, to which we refer 
to as D5j-branes, with 32 D5-branes in each set. (The fact that we need 32 D5-branes in 
each set can be seen from T-duality between these D5-branes and D9-branes.) All of these 
sectors have a well defined world-sheet description in terms of open strings starting and 
ending on D-branes (i.e., they are perturbative from the orientifold viewpoint). On the 
other hand, as pointed out in |16| , the sectors corresponding to the elements flg k and Qrig k 
(k = 1, 2) do not have interpretation in terms of open strings stretched between D-branes. 
In particular, such sectors would correspond to open strings with mixed (that is, neither 
Dirichlet nor Neumann) boundary conditions. The corresponding states cannot be described 
perturbatively. They can be viewed as arising from D-branes wrapping collapsed two-cycles 
in the orbifold [[U|. Such states do not have a world-sheet description, that is, they are 
non-perturbative from the orientifold viewpoint. 

This difficulty is a generic feature in most of the orientifolds of Type IIB compactified 
on toroidal orbifolds. However, there is a (rather limited) class of cases where the would-be 
non-perturbative states are massive (and decouple in the low energy effective field theory) 
if we consider compactifications on blown-up orbifolds ]TjJ. The orbifold we are considering 
here is precisely one of these examples. In fact, we have already pointed out that blowing 
up the Z 3 and (Z 6 )j orbifold singularities is forced on us by the world-sheet consistency 
conditions. Note that the structure of the fixed points in the (Z 6 )j twisted sectors is the 
same as of some of the fixed points in the Z 3 twisted sectors. This implies that blowing up 
the Z 3 singularities automatically results in the corresponding blow-ups in the (Z 6 )j twisted 



sectors. In [11| using Type I-heterotic duality it was shown that all the non-perturbative 



(from the orientifold viewpoint) states decouple in the f2 orientifold of Type IIB on Ai' 



upon blowing up jW' to a smooth Calabi-Yau three-fold. This fact was utilized in (16 
to argue that all the non-perturbative states in the fl orientifold of Type IIB on Ai also 
decouple once we appropriately blow up the Z 3 singularities in Ai. This decoupling can be 
explicitly seen [|ll|] in the heterotic dual (which is perturbative from the heterotic viewpoint 
as there are no D5-branes in the corresponding Type I model) of the fl orientifold of Type 
IIB on Ai' . In particular, on the heterotic side the states corresponding to the flg k and 
flrig k (k = 1, 2) elements of the orientifold group O are charged matter fields arising in 
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the twisted sectors of the heterotic orbifold. These states, which we symbolically denote 
by T a , are perturbative from the heterotic viewpoint. (The heterotic dual is a modular 
invariant orbifold conformal field theory.) The blow-up modes S a of the orbifold also come 
from the twisted sector states in the heterotic model. Moreover, there is a perturbative 
superpotential of the form W = y a bcS a T a Tb + . . . with the Yukawa couplings precisely such 
that the states T a decouple once the blow-up modes S a acquire appropriate vevs (so that 
the orbifold is blown up to a smooth Calabi-Yau three- fold). In [pi it was argued that an 



analogous effect takes place in the Q orientifold of Type IIB on A4. All the naively expected 
non-perturbative (from the orientifold viewpoint) states therefore decouple once we consider 
the corresponding orientifold on a blown-up Ai. 

The above analyses imply that we can use the "naive" tadpole cancellation conditions 
(derived by straightforwardly generalizing those in six dimensional cases of Refs |5|-§] to 
four dimensions) to construct the massless spectrum of the Q orientifold of Type IIB on 
the (blown up) M. orbifold. The tadpole cancellation conditions for examples of this type 
were discussed in fl2|]13|| as a generalization of the cases studied in ||5H8l,|lT||2|l . (Also see 



I^ , p!7lJT5| . ) Here for the sake of brevity we will skip the details and simply state the answer. 

The action of the orbifold group G = {g a \a = 1, • • • , \G\} on the Chan-Paton charges 
carried by the D9- and D5;-branes is described by 32 x 32 matrices 7 ffai g and 7 9ai 5 i , respec- 
tively. These matrices must form a projective representation of G. (A priori they can also 
form a representation of the double cover of G.) Note that the Chan-Paton matrices 719 
and 7i j 5 i corresponding to the identity element of G can be chosen to be 32 x 32 identity 
matrices, and Trf/y^g) = n 9 and Tr(7i )5 J = n 5i , where ng = 32 and = 32 are the numbers 
of D9- and D5j-branes, respectively. (The constraints on the numbers of D-branes arise from 
the untwisted tadpole cancellation conditions.) The twisted tadpole cancellation conditions 
read: 

Tr( 7g , 9 ) = Tr(7 92)9 ) = -4 , (3) 
Tr(7 rij9 ) = Tr(7 grii9 ) = Tr(<y g 2 r . )9 ) = , (4) 

and similarly for the Chan-Paton matrices 7^,5,. A choicef] consistent with the requirement 
that the Chan-Paton matrices form a (projective) representation (of the double cover) of G 
is given by: 

7 9 ,9 = diag(wJi2, uj 2 I 12 , J 8 ) , (5) 
7 rij9 = io-i <g> I 16 , (6) 

where <7j are 2x2 Pauli matrices, and Im is an M x M identity matrix. (The action on the 
D5j Chan-Paton charges is similar.) 

It is not difficult to compute the massless spectrum of this theory. First, consider the 
closed string sectors. Before the orientifold projection we have M = 2 four dimensional 



1 This choice is unique up to equivalent representations. The uniqueness of the choice for TV (7^9) is 
not difficult to see. On the other hand, the uniqueness of the choice for Tr(7 ri)9 ) was argued in 
from the orientifold viewpoint, and was recently shown in [16| from F-theory (l^] considerations. 
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supergravity coupled to h 1 ' 1 = 36 neutral hypermultiplets and h 2 ' 1 = vector multiplets. 
After the orientifold projection we have M = 2 supergravity coupled to h ' l + h 2,1 = 36 [16| 



neutral chiral multiplets. These states are summarized in Table [I] sector-by-sector. Next, 
consider the open string sectors. We have 99, 5;5j, 95j and 5j5j {i ^ j) open string sectors. 
We will confine our attention to the vacuum where the D5j-branes are sitting on top of each 
other at the origin in A4, and take all the 99 sector Wilson lines to be trivial. Then 99 
and 5j5j gauge groups are identical. Thus, the 99 gauge group is SU(6) ® 5p(4) ® U(l). 
(Here we are using the convention where Sp(2M) has rank M .) The massless fields (chiral 
multiplets) charged under the gauge group are summarized in Table |. Note that in each 
SU(Q) subgroup we have 3 chiral generations. (The number of generations in SU(6) is 
determined by the number of chiral representations transforming in 15 of SU(6).) 

Here we should point out that the 99 and 5j5j £7(1) 's are all anomalous ]nj in this 
model. We therefore have the corresponding Fayet-Iliopoulos D-terms. These D-terms 
are cancelled via a generalized Green-Schwarz mechanism ||20|| . There is a subtlety here, 
however. The blow-up modes (that is, the neutral chiral multiplets in the closed string 
sector) transform non-trivially under the anomalous U(l) gauge transformation. (This was 
suggested in and shown in [|ll]] explicitly in the case of the Q orientifold of Type IIB on M.' 
using Type-I heterotic duality.) These blow-up modes (together with the dilaton plus axion 
supermultiplet) cancel the Fayet-Iliopoulos D-terms via acquiring vevs. This is consistent 
with the observation that the orbifold singularities must be blown up [16| as we discussed 



previously. In fact, from Type I- heterotic duality one can see [|TT[] that it is precisely the Z 3 
blow-up modes that are responsible for cancelling the D-terms in accord with our previous 
discussions. As a result, all the £7(1) 's are broken (and must therefore be deleted from the 
massless spectrum in Table |). Moreover, 4 out of the 36 chiral multiplets in the closed 
string sector are eaten in the super-Higgs mechanism. 

This Type I vacuum contains D5-branes and therefore its heterotic dual is non- 
perturbative: under Type I-heterotic duality Type I D5-branes are mapped to heterotic 
NS 5-branes. Let us understand the dynamics of D5-branes in more detail. In particular, 
we would like to study flat directions corresponding to motion of D5-branes. Since we are 
dealing with a four dimensional J\f = 1 vacuum, we need to deduce the (tree-level) super- 
potential in this model. It can be computed using the if-charge conservation rule |TL|. The 
if -charges are given in Table [I]. Here we give the non- vanishing renormalizable terms only 
(we suppress the actual values of the Yukawa couplings): 

W = eij^iXjXk + eijk&iX'jXl + <U^IQ' J Q' J + <^J ri l? J + 

<\\()'Q' + \,r ; ir + p iJ Q jk R k] + 

q'jp'qj . ir j Q'i' j ■ r ,J irii' • ... . (7) 

Here summation over repeated indices is understood. The F-flatness conditions can be read 
off from this superpotential, and together with the D-flatness conditions they determine the 
classical moduli space of the theory. 
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TABLE I. The massless spectrum of the four dimensional Type I Z3 ® Z2 ® Z2 orbifold model 
with A/" = 1 space-time supersymmetry and gauge group [SU(G) ® Sp(4) ® C/(l)] 4 . The number 
of neutral chiral multiplets in each of the closed string sectors is shown. The if-charges in both 
— 1 and —1/2 pictures for states in the open string sector are also given. The gravity, dilaton and 
gauge supermultiplets are not shown. The indices i,j, k, I take values 1, 2, 3. 
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